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I. INTRODUCTION 



Isochronous systems have been quite fascinating ever since Gahleo's discovery of first such 
system - the simple harmonic oscillator. Eminent names such as Newton, Huygens, Jacobi, 
Poincare etc have delved in the investigation of isochronous systems. In the last century 
the focus has mainly been on 2-dimensional cases [4] - as it is related to the Hilbert's 16th 
problem. The other paradigm [1-3] of an isochronous system is the rather unexpected case 
of the potential V{x) — \ {x^ + In a two dimensional dynamical system, isochronicity 
implies that the time period is independent of the amplitude of motion (i.e. independent 
of the initial condition) . An exhaustive study of various isochronous centers of vector fields 
in the plane can be found in [4]. In recent years the seminal work of Calogero and col- 
laborators [5-8] has opened up new and interesting directions in the study of isochronous 
systems. Calogero ct al introduced a simple trick by which one can construct isochronous 
oscillators by modifying equations for ordinary oscillators. And this led to the conclu- 
sion that - "isochronous systems are not rare" [9] .This area of research focuses mainly in 
the construction of isochronous systems and has been widely exploited of late by Ghose 
Choudhury and Guha [10]. In this work however we explore a different angle. Given a two 
dimensional dynamical system we propose to find the conditions for isochronicity by using 
renormalization group technique. 



II. THE METHOD 

The renormalization group approach has been an effective tool for analyzing periodic 
orbits in dynamical systems. First introduced by Chen et al for the purpose of extracting 
asymptotic behavior of the solutions to various differential equations, it has had numerous 
applications over the following years. Recently it has been demonstrated to be an useful 
method for differentiating between center like oscillations and limit cycles for two dimen- 
sional dynamical systems [11, 12]. In fact, that differentiation has been shown to lead to a 
classification of orbits in Lienard systems of different orders [13] i.e. the dynamical system, 
X -\- x"'f{x) + X + g{x) — 0, where n = 1, 2, ■ • ■ . In these Lienard systems, f{x) and g{x) 
are nonlinear functions of x, with g{x) not having any linear term. In this communication 
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we show that the renormahzation group can also be effective for identifying isochronous 
oscillations. 

The methodology of implementing the renormahzation group is as follows: 

1. Identify a fixed point that is either directly a linear center or can be made a linear 
center by suitably choosing a parameter in the problem (this parameter may in certain 
situations have to be introduced). 

2. Set up a perturbation theory around the closed orbit of the center. The orbit will be 
characterized by two constants — the amplitude A and the phase 6, fixed by the two 
initial conditions associated with a second order dynamical system. 

3. The perturbation theory will most likely diverge (due to the presence of secular terms) 
as the separation t — to {t being the current time and to the initial time) becomes large. 
Two renormahzation constants have to be introduced to absorb these divergences from 
the past (to — oo). The renormahzation constants appear in terms of an arbitrary 
time r, which serves to peg the new initial conditions. The amplitude and phase are 
now dependent on r. 

4. The value of x at time t (the current time) cannot depend on where one sets the initial 
conditions and hence — 0. This leads to the flow equations. 



^ = nA) (1) 

If the oscillation is of center variety, then the initial condition sets the amplitude of motion 
and hence dA/dr — 0. For an isochronous orbit of the center variety 

f{A) = 0; g{A) = identically. (3) 

We now present an algorithm for calculating f{A) and g{A). This is best done by considering 
the anharmonic oscillator 

X + u?x = -\x^ (4) 
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The expansion x — xq + Xxi + }?X2 + ■ • ■ leads to the system of equations (different orders 
of A") 



n — Q : xq-\- uj^xq — (5) 
n = 1 : xi+ LiP'xx = —xl (6) 
n — 2 : X2 + ijJ^X2 — — SxqXi (7) 

The initial condition is taken as x(io) = ^ and x(to) = 0. The solution at n = is 

xo = Acosa'(t — io) (8) 

where the initial condition has been absorbed wholly in xq. For n 7^ 0, we need = 
0) = 0. At the order n = 1, 

X\ + uP'xx — —x% — —J^ QQ^ bj{t — to) 

— — — {3 cos uj{t — to) — cos 3u!{t — to)) (9) 

It is cos u!{t — to), the resonating term on the right hand side of Eq. (9) that causes xi{t) 
and hence the perturbation theory, to diverge as t — to ^ 00. The solution for Xi{t) is 

^i(^) = — ~ ^0) sina;(t - to) 

+—-^[cos3oj{t-to)-cosoj{t-to)] (10) 

giving 

x{t) — Acosuj{t — to) — {t — to) smu!{t — to) 

+ [cos3uj{t - to) - cosuj{t - to)] + 0{X^) (11) 

The t — to iii the divergence causing term is now split as t — t + r — to where r is an 
arbitrary time and we introduce two renormalization constants as A = A{T)Zi{T,to) and 
©0 — —i^to — 0(t) + 22{t, to). The constants Zi^2 have the expansions 

= 1 + Aai + A^a2 + . . . (12) 
Z2 = Xh + X%2 + ■ ■ ■ (13) 



and we can write Eq. (11) as 



x{t) = A{t) [1 + Aai + . . . ] cos {ujt + Q + Xbi) 
3XA^ 

— (t-T + T-Qo) sin {Ljt + Q) + ... 

oOJ 

= A{t) cos {uot + e) + XA{T)ai cos {ujt + 6) 

3XA^ 

-A{T)Xbi sin {out + 0) — - r + r - Oq) 

oUJ 

xsm{ujt + e) + .... (14) 

Choosing oi = and bi — —^^-{'t — to) we remove the divergence causing terms (the terms 
containing r — to) and we are left with 

SXA^ 

x{t) = A{t) cos {ut + 0) — {t - t) sin {ut + 0) . (15) 

Since the final solution cannot depend on the arbitrary time scale r, we impose the condition 
(If)* ^ ^' ^' which leads to 

^^0, (xa) 

This is in accordance with what should have happened. The condition, dA/dr — 0, implies 
a center type periodic orbit and dQ/dr = g{A) implies an amplitude dependent time-period. 
We want to draw attention to the fact that dA/dr = is a consequence of the absence of 
smuj{t — to) on the right hand side of the Eq.(9). It will continue to remain zero at 0{X^) 
and above if the right hand side of X2 + uj'^X2, x^ + cu'^x^ and so on never has smuj{t — to) 
term. On the other hand, dQ/dr is determined entirely by the coefficient of cos uj{t — to) on 
the right hand side of the Eq. (9). If this term were absent, then dO/dr — 0, which imphes 
an amplitude independent time-period and hence we would have an isochronous oscillation. 
Thus the signature of isochronous oscillations, dA/dr = and dQ/dr = 0, is produced by 
the absence of any driving term involving sinu{t — to) and cosu){t — to) on the right hand 
side of the successive differential equations for Xn{t). The vanishing of dA/dr and dQ/dr 
gives an isochronous center or else it gives at every order a relation among the constants in 
the problem that has to be satisfied for an isochronous orbit. 
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III. EXAMPLES 
A. A simple case 

As a first check that this statement is vahd, we considered the most well known nontrivial 
isochronous system — the oscillations in the potential V{x) — ^ {^^ + -^)- The equation of 
motion is given by 

x = -x + ^ (18) 

The linear center is clearly at a; = ±1 and it suffices to study the orbit around only one 
of them (there is an infinite barrier at x — and hence the two centers are completely 
decoupled). We shift the origin to the center according to our strategy and write eq. (18), 
in terms of the new variable y — x — 1, 

y^-{l + y) + {l + y)-^ 
= -% + 6y^ - lOy^ + 15?/^ - 21y^ + ... (19) 

For book keeping purposes we introduce a parameter A and rewrite eq. (19) as 

y + Ay^ A6y' - A'lOy^' + A^'lSy^ - A^21y^ + . . . (20) 
Now we can expand y as 

y = yo + Ayi + + X^ys + ■■■ (21) 
and proceed as earher. At diflFerent orders of A we find the following differential equations: 



A°: 


Vo + 42/0 


= 0, 


(22) 


A^: 


yi + %i 


= 62/o> 


(23) 


A^: 


2/2 + 42/2 


= 122/02/1 - lOyl 


(24) 


A=^: 


2/3 + 4^3 


= 12yoy2 + 62/1 - 3O2/02/1 + 152/0 


(25) 



and so on. We work with the initial conditions y = Aq, y = at t = 0. We let 
Xo pick up the initial condition (2/0 = Aq, 2/0 = 0), so that for all subsequent orders — 
Xi{t = 0) = Xiit = 0) = 0. Explicit calculations upto 6th order in A confirms the absence of 
resonating terms (cos fit or sin fit) in the inhomogenous part of the equations at each order, 
estabhshing that upto C(A^), both dA/dr — and dO/dr — and that our condition (3) 
for isochrony is satisfied by this well known center type isochronous oscillator. 
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B. Cherkas System 



We now turn to a well studied and rather interesting example — a dynamical system 
introduced by Cherkas. This system was first studied by Cherkas [18] but the center condi- 
tions found therein were later found to be incomplete. The complete set of necessary and 
sufficient conditions were found in [19] and [20]. The conditions for which the origin is an 
isochronous center was found by [21]. Wc deal with this system in some detail, so that our 
methodology becomes transparent. The equation has the form 

x = y{l + x), (26) 
y — —x — aix^ — a2X^ — a^x"^ — 04(0; + ar,x^)y — aey^. (27) 

We can write this as a second order differential equation, 

X — {1 + x)y + yx, 
— —(1 + x) [x + aix^ + a2X^ + a^x'^ + 04(0; + a^x^)y + Ogy^] 

+ Y^(l-«6), (28) 

which we cast as 

X + X = —(1 + ai)x^ — (a2 + aijx^ — (02 + 03)0;^ — ay,x^ 

-aix{x ^- a^x^) ^- {1 - ao)- , (29) 

= -^x^ - (oi + 02)^^ - (02 + a^)x^ — 03^^ 

—aix{x + 053;^) + ax^ [1 — a; + — + . . . ] . (30) 

where ^ = 1 + oi and a — 1 — a^. The book keeping parameter A that we introduce is 
distributed as follows 

X + X = -X^x^ - X^ (ai + 02) x^ - X^ (02 + 03) x'^ 

—X^a^x^ + a4A,T,i — a4^a^X^xx^ + aXiP' 

x(l -Ax + AV -AV + ...). (31) 

We now expand x as, 

X — xq-\- Xxi + )?X2 + ■ • ■ . (32) 
Using the above expansion in eq. (31), at different orders, we get 



xo + xq^ 0, (33) 

xi + xi — —(3x1 + '^^0 ~ 043^0^0, (34) 
X2 + X2 = —2(3xqXi — (ai + 02)0:0 + 2axiXQ — ax^XQ 

—a^ixiXQ + xqX\) — a^a^XQxl. (35) 

We use the initial conditions x{t = = A and x{t = 0) =0 and impose the condition on 

Xq, so that all subsequent Xn have the initial conditions a;„(0) = Xn{0) — 0. Upto 0{X), the 
solutions are, 



Xq = Aq cos t, (36) 

xi^Al 



a- B a + 13 «4 . ^ , 

H — cos 2t — — sm 2t 



2 6 6 

(37) 



2a — P 04 . 

cos t -\ sm t 



3 3 

Inserting these solutions in the right hand side of eq. (35), we find that the condition for 
the existence of a center (vanishing of the coefficient of the sin t term) is 

a4{/3 - a - as) = (38) 
i.e. 04 = or (3 — a = a^, (39) 

and the condition for isochronicity is 

These conditions hold only for 03 = 0, in the system equation (28). This is because the order 
upto which we have calculated so far does not include the 03 term at all. However, notice 
that the conditions given above are only necessary conditions. They need to be satisfied for 
the existence of isochronous center for — 0, but obviously do not guarantee the existence 
of a center. If we have a center with even 04 = (the case considered by Ghose Choudhury 
and Guha), then for 03 = 0, the isochronicity condition turns out to be 

^a{l + ai) - ^(1 + ai)2 + ^{a, + as) = y - ^, (41) 

where 1 + ai — /3. It is to be noted that all the isochronous cases found with 03 = satisfy 
the condition given in eq.(41). In the interesting special case of a = ^, this condition reduces 



to 



«-2) =4«2- (42) 



In the two isochronous cases [10], we have encountered so far, with a = /3, the resuhs are 
in accordance with eq.(42). The usefulness of such constraints he in the fact that in one's 
search for isochronous systems, one knows in which subspace the search has to be hmited. 
For 03 ^ 0, the calculation needs to be carried through two more orders and at the end of 
a long but straightforward calculation one obtains 

5 A 65 A 355 9 9 65 o 55 n 31 o 
4 48^^ 48 ^ 12 ^ 12 ^ 24 

"12"'^ ^ 48^^ '^ ~ 16'^^"' ^ ""'^ ^ + «2) 

27 1 83 7 

-—a^{ai + 02) - -ap + — + -^(as + ag) 

3 13 
--a{a2 + as) - + «2) + ^(oi + aa)^ 

-^«3 - ^ = (43) 

Once again all the 03 7^ cases noted by Ghose Choudhury and Guha satisfy the above 
constraint. 

When 03 = and a — P the above constraint takes on a much simpler form given by, 

1 2 / \ 19 2 1 1 / \ 

a (oi + aolH a H — aao aiai + ao) 

48 16 ^ ^ 192 4 16 ^ ^ 

+ ^{ai + a,f - I = 0. (44) 

For Cherkas system when 03 = and « = /3 all the isochronous systems must satisfy the 
above condition. Keeping in mind the fact that a = /3 = 1 + ai, the above constraint gives a 
relation between oi and 02 much like the constraint obtained at previous order, given by the 
Eq. (42). So we have two equations for two parameters which must be satisfied for isochrony 
and thus the two constraints ( Eqs. (42) and (44)) constitute both necessary and sufficient 
conditions for an isochronous Cherkas system, with 03 = and a — 13. Now eliminating 02 
between Eqs. (42) and (44) we are left with a cubic equation in a, given by 

Aa^ - 24q;^ + 45q; - 27 = 0. (45) 

There are only two roots of the above cubic equation i.e. a = 3 and a = 3/2 (which is a 
repeated root). These two cases corresponding to the parameter values of: a = 3, 02 = 1 



and q; = 3/2, 02 = are the only possible values for which Cherkas system (when 03 = 
and q; = ^) is isochronous. 

To further establish the efficacy of our method we consider a set of cubic systems as defined 
by Bardet et al given by, 



X = —y{l + aix) + a2X^ + 032;^, 

ij = X + bsx^ + b^x^ — y{b2X — b^x^) — biy^. 

The necessary condition for the isochronous centers as found by us is 

(2a2 - 62) 



ai + bi 63 + ai 3 , bi ai 

—3 ^ + 4^^- 4+y 



and 



+02(261 + oi) — 64 — 803 = 



^^i±M(ft3 + + 200(63 + 61 + 2o0 - + 
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x{7bi - 563 + 2oi) - ^(^5 + aib3 - 0262) 

_Oi(oi + M +1^202-62)^ = 
4 12^ ^ 



(46) 
(47) 



(48) 



(49) 



All the system mentioned in Theorem 4.1 of Bardet et al satisfy the criterion given above. 



C. Ricatti Equation 

We finally turn to a variant of the Ricatti system of Type-II, given by 

X + kxx + -^x^ + axi — 0. (50) 
9 

Oscillators having equation of motion of the form x + f{x)x + g{x) — are called Lienard 
type oscillators. For our example f{x) — kx and g{x) — {k'^/9)x^ + axi. This oscillator has 
been widely treated in the literature and has been referred to as with various names. The 
oscillator may be interpreted as a cubic anharmonic oscillator with a damping type nonlinear 
force kxx. The equation has otherwise been called the generalized Emden-type equation [24] 
and has been widely studied in the literature over the last two decades. The equation arises 
in a variety of physical problems — like the modeling fusion of pellets [25] , one dimensional 
analogue of Yang- Mill's boson gauge theory [26] , equihbrium configurations of a spherical 
cloud subject to the laws of thermodynamics etc. [27]. Recently Chandrasekhar et al [23] 
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have shown that for o; > 0, this oscillator admits non- isolated periodic orbits with the 
unusual property that frequency of oscillation remains independent of amplitude and same 
as that of the linear oscillator. For our purposes we will show that the isochronous periodic 

solution that the eq. (50) admits, also follows the criterion we have put forward in this 
paper, i.e. dA/dr = dQ/dr = 0. Assuming that k is small we expand x in powers of k as 
follows: 

X — Xo + kxi + k'^X2 + k^X3 . . . (51) 

We work with the initial condition as a;(to) = A and i;(to) = and absorb the initial condi- 
tion wholly at zeroth order as we had done earlier. Proceeding henceforth according to our 
prescription we calculate the flow equations perturbatively. We find, as per our expectations 
the flow equations turn out to be dA/dr — dQ/dr — 0, right upto order 0{k^) and we have 
an indication of an isochronous type oscillator. Numerical simulations of the oscillator 
represented by eq. (50), reveal that for values of a{= Vl^) > 1, one obtains isochronous 
non-isolated orbits and so is the conclusion of the work by Chandrasekhar et al [23]. Using 
our RG approach we have shown here that for a oscillator to execute isochronous oscillations 
it must obey the criterion we have proposed, i.e. both dA/dr and dQ/dr must be zero. 



IV. CONCLUSION 

The above example is a rather interesting case where even after addition of nonlinearity 
the oscillator frequency remains independent of amplitude and same as that of the lin- 
ear harmonic oscillator. This opens up the possibility of identifying more such nonlinear 
oscillators which admit oscillations with amplitude-independent frequency. One can look 
to identify suitable nonlinearities which, added to the linear harmonic oscillator case, do 
not generate resonating terms at different orders of perturbation. Thus our method can 
in principle provide the necessary conditions for finding oscillators having this property. 
Moreover we can ask the question what is the most general case of a nonlinear oscillator 
of the form x — —x + f{x,x) that can admit solutions exhibiting isochronous center-type 
oscillations. We believe our method has important implications in developing nonlinear 
systems exhibiting isochronous oscillations. The obvious limitation of our methodology, is 
being perturbative, it provides only the necessary and not necessarily sufficient conditions 
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for isochronicity. However due to its simplicity it allows us to get an immediate result for 
whether a particular set of nonlinearities will lead to isochronicity or not and thus also indi- 
cates the tinkering necessary to get an isochronous system. Further this method can serve 
as a robust way to perturbatively construct the so called period functions which must vanish 
for isochrony. Much of the work in the field of isochronous oscillations involves finding a 
Hamiltonian description of the system and going on to find conditions for integrability or 
super-integrability. However the technical difficulties involved have prevented a complete 
solution (e.g. for the Cherkas system) at one shot. Besides many cases can't be immedi- 
ately reducible to Hamiltonian systems and in most cases it proves rather difficult to find 
first integrals when nonlinearities arc involved. Our method however is not limited by such 
considerations and can be used effectively as a first probe for isochronous dynamical systems. 
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